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Discussion of 
“THE DEVELOPMENT OF THE TURBULENT BOUNDARY LAYER 
ON STEEP SLOPES” 


by William J. Bauer 
(Proc. Sep. 281) 


G. HaLBronn.'*—A large number of experiments have been performed by 
Mr. Bauer which were concerned with boundary-layer development on spill- 
ways having various conditions of slope and roughness. This phenomenon was 
previously studied, both theoretically and experimentally.'® It is of interest 
to compare the author’s results with those obtained from a study of Eq. 3 
when u is assumed to follow one of the relationships valid for the velocity dis- 
tribution in the case of uniform flow within a pipe. The mathematical deri- 

vation can be evolved when: 


Total Head ° 
SS (a) The region under study is 
located away from the spillway 


crest. In this case 


U= V2g(d + Ho — yocos a). (23a) 


can be expressed approximately as 
U = ¥2g(d + Ap)... (23b) 


The terms in Eqs. 23 are defined 
in Fig. 12. 

(b) The spillway, away from 

Fig. 12.—Dertnrrion SKETCH the crest, has a constant slope a. 

It can be seen that the boundary- 

layer thickness is independent of the flow near the crest and is dependent only 

on the local conditions. Therefore, the results are general. 

When the spillway is smooth, Eq. 4 can be used if 


Letting 
56 
[To 
p 
Eq. 4 can be written as 
6 0.278 
R, z (26 
in which 


Highway and Bridge Engr., 6th Elee. Dept., Grenoble, France 


» “Etude de la mise en régime des écoulements sur les ouvrages Aforte pente. Application au probléme 
de lentrainement d'air,"” by G. Halbronn, Houtlle Blanche, 1952, No. 1, pp. 21-40; No. 3, pp. 347-371; 
No. 5, pp. 702-722. 
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Therefore, from Eq. 3, 
R, = 4.35 e (22 — 2.82 + 4.32) 


The variation of 6/z, as a function of R,, is given in parametric form by 
Eqs. 26 and 28. This variation can be approximated by 


= 0.039 


for 8.5 < z < 12.5 (0.00205 < c, < 0.00443). If it is assumed"* that 


then 
6 = 0.158 


The curves corresponding to Eqs. 29 and 31 are shown in Fig. 13. For 
small values of the Reynolds number, the curves are in close agreement. This 


Symbol Boundary Slope 


Celluloid, Haibronn 40° 
Concrete, Hickox 55° 
Glass 
Glass 
Glass 
Screen 
Screen 


Glass 
Glass 
Glass 


Screen 
Screen 


Screen 


Value of the Ratio 


Equetion 31 Equation 29 


10° 
Value of the Reynolds Number 


4 
Fia. 13.—Tue Ratio FUNCTION oF For SmootuH BouNnDARIESs 


agreement indicates that the application of a logarithmic law to the experi- 
ments'* is proper. In the range of Mr. Bauer’s smooth-boundary experiments, 
Eq. 31 accounts for the observed increase of the thickness of the boundary 
layer. This agreement is expected as velocity profiles are better approximated 


Meters Developments in Fluid Dynamics," by 8. Goldstein, Clarendon Press, Oxford, England, 
52, p. 339. 
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by an expression such as 


than one in the form: 


although, in Fig. 2, 7 would be equal to 1/7.5. It was also observed that, for 
large values of the Reynolds number, the velocity profiles do not correspond 
to the logarithmic law. Because these larger Reynolds numbers occurred 
downstream from the critical point, it was assumed!* that the poor agreement 
was caused by a physical change in the flow. The author’s experiments show 
that this is not true. 


Us 


It should be noted from Fig. 13 that, when the Reynolds number 3 


exceeds 10°, the smooth-boundary experimental points depart from Eq. 31 and 


0.12 


34, a = 20° 


0.06 


Boundary Layer Thickness &, in Feet 


Glass 20° 
0.04 ° Glass 40° 
Wi ° Glass 60° 
ve Screen 20° 
wr 6 Screen 40° 
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Glass 20° 
Glass 40° 
Glass “60° 
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10 12 


Distance z, in Feet 


Fic. 14.—Tuickness or THE BounpARY LAYER AS A FUNCTION OF z 


more closely agree with Eq. 29. This value (10°) is the upper limit for the 
validity of Blasius’ law. 
Eq. 31 can also be written as 


0.0106 


6 (ft) = yr 


(34) 


In Fig. 14 the curves corresponding to Eq. 34 for a = 20° and a = 60° are 
shown. As a result of the exponent 0.1 in Eq. 34, 6 is approximately inde- 
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pendent of 8, but this fact has no special physical significance. Curves repre- 
senting Eq. 29 in the form: 
0.0104 


6 (ft) == z 
are also shown in Fig. 14. These curves are in superposition. 

In Fig. 15 there is shown the relationship between c, and R;. This rela- 
tionship is obtained by letting u = U and y = 6 in Eq. 30 to result in 


—\3 
To 


= 0.0463 (36) 
The smooth-boundary experimental points in Fig. 15 are farther from the curve 
of Eq. 36 than the same points are from the curve of Eq. 16 which corresponds 


Symbol Boundary 
Glass 


$6 8 10 
Value of R, x10~* Value of the Reynolds Number x 10 ~* 


Fig. 15.—Tue Revationsuie Between Fic. 16.—Tue SHape PaRAMeTER A8 A FUNCTION 
cy AND Rg or THE REYNOLDS NUMBER FOR 
Smootu BounDARIES 


to the logarithmic law. It is difficult to understand the reason for this devia- 
tion because the situation is reversed in Fig. 13, and both figures are only two 
different methods of representing the same result, obtained by substituting 
either Blasius’ law or the logarithmic law in the von Ka4rm4n equation. 

It is also of interest to study the values of the shape parameter A. The 
logarithmic law leads to 

= 1.765 76 

and 


bn = 1.765 Ve, (1 — 3.54 Vey) 


Then, from Eq. 16, 
A 


1.15 logie AR + 0.645 


Blasius’ law leads to 6, = 5, = 7 6, and A = 2. The corresponding 


curves are plotted in Fig. 16. The experimental points are somewhat scat- 
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Screen 40° 
seal 


tered; nevertheless, the value A = 1.285 appears to be valid for large values 
of R. The curves for Eqs. 18 and 38 follow the general trend of the experi- 
mental points more accurately, although with a systematical displacement. 
It is worth while to note that the curve for Eq. 38 and the curve for the empiri- 
cal law (Eq. 18) are in close agreement. 

In Fig. 16 the experimental points shown by Mr. Bauer in Fig. 10(a) have 
been plotted exactly as shown in Fig. 10(a) with the value of the abscissa 
changed to “Value of the Reynolds Number X 10~*.”’ Also, in Fig. 16 the 
curve for Eq. 18 is shown on the basis of the “new” abscissa. This value 
for the abscissa is correct, rather than Mr. Bauer’s ‘‘Value of the Reynolds 
Number X 107°,” as can be seen by introducing into Eq. 18 the value of R 
equal to 10,000 so as to result in A equal to 1.287. It is true that introducing 
a value of R equal to 10,000,000 into Eq. 18 results in A equal to 1.380, but 
Messrs. Tetervin and Lin established Eq. 18 only for values of R_ between 
1,500 and 40,000. 

When (k/v) ¥70/p is greater than 5, dimensional analysis shows that 
u/W To/p depends simultaneously on (y/v) V and y/k. For (k/v) t0/p 
greater than 71, and for the special type of roughness used by Mr. Nikuradse, 
it has been shown‘ that for uniform flow in a pipe Eq. 5 is valid for the velocity 
distribution. Therefore, 


Applying Eq. 3 to Eq. 39 results in 


e (: 2.24 ) (40) 


Eqs. 39 and 40 give the variation of 6/z as a function of r/k (shown in Fig. 17) 
which can be approximated by 


for 6 <z < 9 (0.00395 < c; < 0.0089) or, in a dimensional form: 


The curves for Eq. 42 are shown in Fig. 14 for k = 0.009 ft, k = 0.0021 ft, 
and k = 0.00053 ft. This last value of k is seen to be in good agreement with 
the experimental data. The considerable variation of the values for k given 
in Fig. 14, in Fig. 9, or by the tests involving uniform flow (although, for a 
screen made with wire of an 0.0008-ft diameter, an equivalent sand roughness 
of 0.009 ft seems unexpected) leads one to believe that Eq. 5 does not agree 
with the observed velocity-profile distribution. In fact, if cy is determined 
from the experimental points of Fig. 3 and from Eq. 3 and is then used to plot 
u/~ 7o/p as a function of y, a single curve does not result, as it would if Eq. 5 
were applicable. It would be most interesting to investigate Mr. Bauer’s data 
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on rough bottoms by the method described by Boris A. Bakhmeteff.'"’ By 
plotting [(u/V70/p) — 5.75 logio (y/k)] as a function of log [(k/v) 
it may be possible to determine the velocity distribution corresponding to the 
observations. Examination of Figs. 8 and 9 appears to show that ¥ 71o/p 
depends on both y/k and (y/v) To/p. 

Another manner of solving the problema might be by the use of prototype data 
obtained by the simple method followed by Mr. Hickox.'* Such results (in- 
cluding a measured value for 5) have been published'* by V. Michels, A.M. 
ASCE, and M. Lovely. These results, obtained at the Glennmaggie Dam in 


Australia, are as follows: 


4, in feet 4, in feet 
z, in feet (measured ) (computed) 


0.127 
0.322 
0.525 


Plotted in Fig. 17 (with k = 0.005 ft for concrete) these points fall between 
the author’s design curve and the curve of Eq. 41. The foregoing listing also 
shows the values of 6 computed by the method used by Mr. Bauer in applying 


Symbol Boundary Slope Cu Ft per Sec 

Concrete, Michels per Ft 
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Screen 20° 
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Fie. 17.—Tue Ratio Function or THe Revative Rovenness ror Roven Bounparies 


Mr. Hickox’s data. However, it is noticeable that the measured values are 
larger than the computed ones. This is because the measured values account 
for the corrugation of the surface at the critical point. 

It is also true that the Manning coefficient for concrete can be expected to be 
between 0.013 and 0.016, corresponding approximately to 0.004 ft <k <0.016 ft. 
A change in the value of k which is used to plot the points in Fig. 17 results 

The Mechanics of Turbulent Flow,” by Boris A. Bakhmeteff, Princeton University Press, Prince- 
ton, N. J., 1936, p. 90. 


“Some Observations of Air Entrained Flow,” by V. Michels ant Lovely, Proceedings, 
Minnesota International Hydraulics Convention, Minneapolis, Minn., 1953, p. 4 
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in a constant translation parallel to the abscissa. For these reasons, it can 
be concluded that, whereas the general trend of the design curve is known, 
its position is not accurately determined. 

From Eq. 34 it can be seen that 


dé 0.0074 


For z equal to 10 ft, d5/dz equals 0.0037. The divergence of the boundary 
layer is then 0° 13’. In the same manner, from Eq. 42, 


0.154 


For k/z equal to 0.0005, the divergence is 0° 40’. These low values for the 
divergence help one to understand why velocity distributions, valid in uniform 
flow, can also be applied to this case. 

Mr. Bauer presents the experimental data on a physically sound basis. 
The data allow comparison with the results of a theoretical method of approach 
which was previously developed.'* It is therefore possible to conclude that 
this method is accurate enough to predict the development of the boundary 
layer on smooth spillways. On rough spillways, it remains a difficult problem 
to decide which value of roughness to select. This problem, however, is not 
confined merely to this case; it will be solved when the mechanism of turbulent 
flow on rough boundaries is clearly understood. Fortunately, as pointed out 
by the author and as shown in Fig. 14, a large variation of k has little influence 
on the boundary-layer thickness. This circumstance enables the solution of 
design problems by the method described by Mr. Bauer. 


Wituiam J. Baver,” J. M. ASCE.—An able analytical appraisal of the 
writer’s data has been made by Mr. Halbronn. Based on the needs of the 
researcher rather than the practicing engineer or designer, his discussion is 
an excellent example of how valuable sidelights can be derived from the lab- 
oratory data. In general, Mr. Halbronn’s findings support the conclusions. 

Mr. Halbronn has developed two expressions (Eqs. 29 and 31) for 6/z as a 
function of the Reynolds number R,. Essentially, Eq. 29 is based on the 
assumption of velocity distribution of the form of Eq. 33, and Eq. 31 is based 
on one of the forms given in Eq. 32. Testing the hypotheses thus developed 
in the light of the data shown in Fig. 13, Mr. Halbronn concludes that an ex- 
pression similar to Eq. 32 better approximates the shape of the velocity profile. 
It is felt, however, that the actual shape of the velocity profile within a turbu- 
lent boundary layer cannot be accurately approximated by simple mathe- 
matical expressions of the types considered (Eq. 32 and Eq. 33), except over 
specified limited ranges of the Reynolds number. It should be noted that in 
Figs. 2 and 3 there is a consistent deviation (however slight) of the experi- 
mental points from the straight line which represents an expression of the form 
of Eq. 32. The deviation increases with increased values of the Reynolds 


' Senior Designer, Harza Eng. Co., Chicago, II. 
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number, demonstrating the inadequacy (in this range of the Reynolds number) 
of the simple mathematical approximation. 

Mr. Halbronn notes that in Fig. 2 the value of the exponent j of Eq. 32 
would be 1/7.5. Actually, 7 in Fig. 2 decreases gradually as the Reynolds 
number increases, varying between 1/7.3 and 1/8. Similar trends were ob- 
served in the remainder of the data which were not published. Mr. Schlichting 
states® that beyond the range of validity of the Blasius law (the 1/7-power 
expression) ‘“* * * a better approximation is obtained by the power 1/8, 1/9, 
or 1/10 instead of 1/7." The experimental data seem to confirm Mr. Sch- 
lichting’s observation. 

In analyzing the data for c,, Mr. Halbronn has assumed a 1/7-power velocity 
distribution and has derived Eq. 36, an expression for c; in terms of Rs. On 
the basis of the previous comments, it is of interest to repeat the procedure 
for an 1/8-power velocity distribution. This procedure results in 


cy = 0.0427 


The curve for Eq. 45 lies about as far above the data in Fig. 15 as Eq. 36 lies 
below, showing that the best fit would be associated with a value for j of be- 
tween 1/7 and 1/8. This fit agrees with the previous indications as the ob- 
served values of j are between 1/7.3 and 1/8. It can be concluded, therefore, 
that the lack of agreement between Eq. 36 and the data is not caused by the 
inappropriateness of the form of the equation which has been used to approxi- 
mate the velocity distribution, but rather to the use of an inappropriate value 
of j. 

The excellent agreement between the data and Eq. 16 is surprising, as Mr. 
Halbronn has indicated, as this equation is based on the apparently inappro- 
priate logarithmic velocity distribution. It is believed that this agreement is 
merely coincidental. The apparent agreement between the data and Eq. 31 
in Fig. 13 is also somewhat coincidental as Eq. 31 is based on a slightly erroneous 
value of j, whereas the definition of 6 in Fig. 13 does not agree exactly with 
the definition of 6 used to derive Eq. 31. These two approximations have 
apparently tended to compensate for each other. 

A misunderstanding has arisen concerning the significance of Fig. 10(a). 
Mr. Halbronn has suggested that the error lies in a mislabeling of the loga- 
rithmic abscissa scale. Actually, the proper abscissa scale for Fig. 10(a) 
ranges between 3 X 10° and 5 X 10’ and is the value of R, = U z/v and not 
R;, = U6,/v. What the writer failed to explain was that Eq. 16 was first 
modified by a substitution of 


and then plotted in Fig. 10(a). Eq. 46 was determined from experimental 
data” (not published herewith). 

In order to analyze the development on the rough bottom, Mr. Halbronn 
used the von Karm4n-Prandtl velocity-distribution equations for uniform flow 


~The Development of the Turbulent Boundary Layer on Steep Slopes,"’ by William J. Bauer, dis- 
sertation presented to the State Univ. of lowa, at Iowa City, lowa, in August, 1951, in partial fulfilment of 
the requirements for the degree of Doctor of Philosophy. 
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in pipe. Using Eq. 3 and Eq. 5 for the rough pipe (and making several nec- 
essary assumptions), he derived an expression for 6 as a function of z and k, 
plotting it in Fig. 14 for three fixed values of k. Inasmuch as the assumed 
logarithmic velocity distribution is even further from the truth in the case of 
rough boundaries than in the case of smooth boundaries, the subsequent analy- 
tical results are of questionable value. In Fig. 3 are illustrated velocity pro- 
files approximated by equations of the type of Eq. 32 wherein j ranges between 
1/4.23 and 1/4.86. These values are far from the 1/7-exponent and 1/8- 
exponent which give better (yet inadequate) agreement with the logarithmic 
distribution. 

The writer is pleased with the insertion of additional prototype data into 
Fig. 17, as well as with the tabulated comparison between measured and com- 
puted values of 6 on the Glennmaggie Dam. It is agreed that the position of 
the design curve in Fig. 17 is not accurately determined (k remains in doubt). 
It should be noted (as Mr. Halbronn noted) that the range of doubt is 
sufficiently small to permit solution of design problems by the method illus- 
trated in Fig. 11. Conservatism would merely indicate the assumption of a 
reasonably large value for k. 
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Discussion of 
“AN EXPERIMENTAL INVESTIGATION AND LIMIT ANALYSIS 
OF NET AREA IN TENSION” 


By W. G. Brady and D. C. Drucker 
(Proc. Sep. 296) 


W. H. MUNSE,! A.M. ASCE.—It is gratifying to see a new approach applied 
to a problem which the engineering profession has been studying for many 
years—that of determining the net area of a tension member. The writer, 
however, hopes that the analysis and tests presented in the paper by Messrs. 
Brady and Drucker will not be considered as the answer but only as a part of 
the answer or a starting point upon which it may be possible some day to build 
a better answer than is now available. These tests, as well as many of the 
previous tests on net section, have been conducted on flat-plate specimens 

and then analyzed in terms of design rules which are used for riveted or bolted 
connections in which the load transfer may differ greatly from that of the flat 
plate. The engineers need an accurate analysis or analyses for the behavior 
of tension members but will want assurance that the theory is applicable to 
realistic members. 

The limit load theorems presented in this paper provide a new viewpoint 
insofar as the behavior of tension members is concerned. Consequently, it 
would have been helpful if the tests had been carried to failure and the rela- 
tionships between the limit load and ultimate load-carrying capacity of the 
plates determined. This information would have made it possible to compare, 
at least in part, the result of these tests with the many tests reported in the 
literature. 

In some structures, general yielding across an entire member may produce 
failure of the member or, at least, exceed its limit of serviceability; in other 
instances initial yielding or possibly, ultimate load-carrying capacity may be 
the limiting factor of serviceability and a more reasonable value upon which 
to base the working load. Thus, if the material is to be used efficiently, the 
mode of failure must be considered in any determination of working stresses 
and not only the limit load. 

Consider, for example, a long tension member in a bridge. When general 
yielding occurs in the connections of such a member, this yielding is usually 
restricted to the first row of fasteners in the connections at either end of the 
member: a relatively small portion of the member. The effect of this plastic 
deformation in the connections on the overall deformation of the entire mem- 
ber may be small and have little effect upon the usefulness of the member; in 
all probability, it would be a general yielding of the entire member which 
would be the limiting factor for such a structure rather than the yielding at 
the connections. However, even after the entire member has been subjected 
to a general yielding, the member would still retain a considerable amount of 
additional load-carrying capacity. 

In recent tests conducted at the University of Dlinois on 18 in. I-sections 


1. Research Associate, Prof. of Civ. Eng., Univ. of Dlinois, Urbana, Ill. 
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with 1/2-in. gusset plates, it was found that the distribution of load in the 
members was greatly different from that generally obtained in flat plates. 
This results from the fact that the load is taken out of the member through 
only the flanges while a large percentage of the area of the members is con- 
centrated in the web. However, strain measurements on the webs of the mem- 
bers provided data which can be correlated with the limit analysis study by 
the authors. A general indication of the behavior of the members may be ob- 
tained from the following tabulation. 


Loading Conditions Average Gross Stress, psi 


Initial Yielding 15,000 
General Yielding in Joint (Critical Section) 36,400 
General Yielding in Gross Section 37,600 
Ultimate Strength of Member 50,200 


From the test data, it was found that general yielding in the joint (when the 
yield point strain was attained in the center of the web) occurred when the 
stress ratio T*/T, or yield-strength efficiency for the member was 97 per- 
cent. However, the members failed at a test efficiency of approximately 74.5 
percent. It appears that parts of the section may have been strained into the 
work-hardening range when the center of the web at the critical section 
reached the yield point strain. These data suggest that several analyses may 
be necessary to explain completely the behavior of tension members: an elas- 
tic analysis, a limit analysis, and an empirical ultimate strength relationship. 

It is realized that the above discussion raises questions concerning the ap- 
plicability of limit analysis totension members. Nevertheless, the writer hopes 
that the authors will find it possible to pursue their study further and to ex- 
pand it to include actual joints of various types as well as such other factors 
as variations in, pitch, gage, material properties, and the preparation of the 
fastener holes. By determining, in such tests, the limit loads as well as ul- 
timate load-carrying capacity it may be possible to evaluate this information 
in terms of the limit analysis theorems as well as other concepts based on 
the ultimate strength or elastic properties of the materials. 


P. P. BIJLAARD,2 M. ASCE.—Calculations similar to those made by the 
authors, were published in Dutch by the writer in 1931, using the maximum 
shear stress yield criterion,? as well as the Mises criterion.4 These calcu- 
lations covered the subject more in detail than the authors did. A summary 
in English was given in a footnote on page 35 of reference 5. In addition to 
the calculations the writer made elaborate experiments.3 Since both refer- 
ences 3 and 4 are in Dutch, the principal results will be given below. It will 
be shown in which situations the authors obtained essentially the same results 
as those of the writer’s earlier investigations, and in which cases the two 
differ. 


2. Prof. of Structural Eng., Cornell Univ., Ithaca, N. Y. 

3. “Resistance of a section of a plate between rivets in an oblique row, calcu- 
lated from the Mohr criterion of plasticity” (in Dutch), by P. P. Bijlaard, 
de Ingenieur (the Hague, Holland), Vol. 46, 1931, pp. B.47-B.57. 
“Resistance of a weakened oblique section of a plate subjected to tension, 
calculated by the plasticity criterion of Huber-Hencky” (in Dutch), by P. P. 
Bijlaard, de Ingenieur (the Hague, Holland), Vol. 46, 1931, pp. B.249-B.251. 
“Theory of local plastic deformations”, by P. P. Bijlaard, Publ. Interna- 


tional Association for Bridge and Structural Engineering, Vol. 6, 1940-41, 
pp. 27-44. 
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The writer showed that minor plastic deformations in a weakened section 
(AB in Fig. 1 of ABB! A! in Fig. 3) result in a state of stress which produces 
maximum resistance against deformation before the section fails. This as- 
sumes an elasto-plastic stress-strain diagram of the mild steel type. This 
situation applies if the maximum shear stress yield criterion is used and, at 
the same time, plastic deformation is assumed to occur by sliding along the 
planes of maximum shear stress.° It also applies if the Mises criterion is 
used and, at the same time, plastic deformation is assumed to occur quasi- 
isotropically, that is, in the same way as the elastic deformation, but with a 
varying plastic modulus and with a Poisson’s ratio of one half.49 

If the maximum shear stress criterion is used for an oblique section like 


that of Fig. 1, one obtains a limit curve for the total stress yu = ag” +T 
which is given by Fig. 2 for any angle a. For @ = 0 the total stress is U *9 
= ys * 2k, where Oys is the yield stress in tension and k is the yield stress 
in pure shear. The limit curve consists of a quarter circle with its center at 
o =k for a ¢ 45° and a straight line 7 = k for a 245°, Fora >45° the state 
of stress in the oblique section AB at the limit load is given by a circle that 
touches the line PN at the end point Q of the pertinent vector “ap which rep- 
resents the total stress uw on the weakened section AB (Fig. 2). On the other 
hand, the pure tensile stress assumed by the authors in obtaining their Equa- 
tion (15) for ao 245° would lead to a limit curve given by the dashed quarter 
circle NO which results in much lower limit loads than the limit curve PN. 

At the same time? the writer showed that, if the plate is subjected to ten- 
sile forces P, subsequent to fully plastic deformation in the section AB, the 
relative movement of the elastic part of the plate below AB with respect to 
that above AB is given by the dashed arrow D, denoted as the direction of 
yielding. For a 445° as in Fig. 1, where the planes of maximum shear stress 
are not perpendicular to the middle plane of the plate, the arrow D subtends 
an angle w*a with the axis of the plate. If a 245° the planes of maximum 
shear stress are perpendicular to the middle plane of the plate, so that the 
direction of yielding coincides with AB, and w= 90° - a, 

In a symmetrical case, as in Fig. 3, it is not necessary that the forces F, 
transferred by the oblique sections AB and A1pl, are in the direction of the 
total force P, acting in the axis of the plate. If the forces F deviate by an 
angle y from the axis of the plate, they will subtend angles @- y with the 
normals to the sections AB and Alp}, Then, from Fig. 4, the total stress yw 
in AB and A1B! is given by the vector OR and the effective stress in the 
direction of the force P by “Me = uw cos y (Fig. 3), given in Fig. 4 by the vector 
OS. It follows directly from Fig. 4 that “~~ * OS will be as large as possible 
if the perpendicular RS to OS is tangent to the limit curve MNP. In that case 
CR is parallel to OS and the angle MOR = 2 - Y *q@/2, so that also y *a/2. 
Hence the limit curve for the effective stress “,. for symmetrical sections is 
given by the locus MSP of the feet S of the perpendiculars RS, which was de- 
noted as the pedal curve to the curve MNP. From Fig. 4M, * k (1 +cosa@) 
(reference 2, page B.54), so that the maximum force which a section AB or 
Alp! in Fig. 3 can take is 


P/2=ktL (1+ cosa) (1) 


where L ® length, t * thickness. 
The author’s Equation (31) is identical with this equation, since their angle a@ 
is the complement of the writer’s angle a. 

Comparison of Figs. 1 and 3 shows why the forces F, acting on each of the 
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sections AB and Alp! will establish themselves in such a way that y * a /2. 
In that case the forces F subtend angles a - y * a /2 with the normals to the 
sections AB and A!B!, Hence from Fig. 1 the directions of yielding D will 
subtend angles w* a /2 with the forces F, so that from Fig. 3 the vectors D 
are parallel to the axis of the plate. It is obvious that the directions D must 
be parallel to each other, since the elastic lower part of the plate in Fig. 3 
moves in one direction only. Further, from symmetry, D must be parallel to 
the axis of the plate. 

The case of the plate between oblique rows of holes is somewhat more in- 
tricate, since one has to determine which section between the holes will fail 
at the smallest limit load F (Fig. 5). In Fig. 2, the limiting total stress yu 
for a single oblique section is bounded by the curve MNP. The oblique sec- 
tions which can take the same limit load as a normal section AA! (normal to 
the force F), with cross section A, must have cross sections 2k A/u . Using 
u“ from Fig. 2 it follows that all those equivalent oblique sections through a 
point (Fig. 5) are bounded by the curve ANA!N!, consisting of 2 straight lines 
of lengths equal to AA! and two semi circles with diameters equal to AA!. 
For example, any oblique section LL, subtending an angle o@ with the section 
AA}, which is normal to the force F, has the same limit load as the section 
Aas}, Consequently that section LL between the holes which can take the 
smallest load is obtained by placing the limit curve ANAIN!, using an appro- 
priate scale, between the holes, so that it touches the perimeters of the holes. 
Then the section connecting the two tangent points, like LL in Fig. 5, is the 
weakest section. It was shown in reference 2 that the line LL must be 
straight. Hence in Fig. 5 the limit load for the section between the holes is 
equal to the limit load for the normal section AA}, or equal to 2 k (AA!) t, 
where t is the plate thickness. 

For small values of § the holes will touch the straight parts of the limit 
curve ANA! NI, so that then, with the present notations (Equation (1) of ref- 
erence 1), 


AA! = - d (2) 


For larger angles 8 the holes touch the circular parts of the limit curve 
ANA!N , 80 that (Equation (2) of reference 1) 


Q? 4 2 _ 4? 
Assuming, for purposes of design, the load F to be uniformly distributed over 


the section MM, in the line connecting the centers of the holes, the total limit 
stress my (in the direction of F) on this section would be 


The limit curve for the oblique limit stress 4 my is given graphically in 
Fig. 6 for the most unfavorable case when \/d is as small as possible, which 
was assumed to be Q/a = 3. This curve was denoted as the upper limit curve, 
since it was obtained, as shown above, by assuming the stress distribution in 
all planes to be as favorable as possible. Obviously these most favorable 
stress distributions differ for adjacent sections, so that they cannot always 
occur simultaneously. 

In order to find a lower limit curve for several angles § (Fig. 5) a body of 
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equal plastic resistance was drawn between the holes, also assuming N/a = 3. 
This body was given such a form that practically everywhere within that body 
the state of stress, @;, 9, was at the yield stress, and such that the resist- 
ance of its effective section AB to the forces F was as large as possible (Fig. 
7). The stress is uniformly distributed over this section AB, which subtends 
an angle a =8 + § with the normal to the line of action of the forces F. In 
Fig. 7 the angle 8 + 6 is chosen as tan~! 2.5, so that from Fig. 2, where 
Hap is the total limit stress on AB, the principal stresses are 01 = 0.7 Cys: 
02" -0.3 dys. Fig. 7 speaks for itself. The load F that can be transferred 
between the holes is F = (AB) tap, So that the stress uj, assumed to be 
uniformly distributed over the section MM = ({ - d) t (Fig. 5), is (Equation 
(7) of reference 1) 


* [(AB)/(MM)] Hap 322 (5) 


Since a more favorable stress distribution than that of Fig. 7 may be possible 
and since the material outside the body of equal resistance will restrain the 
plastic deformation, the stresses ,,¢determine a lower bound for the stress 
Um- By considering sections between holes for various angles § the limit 
curve in Fig. 6, denoted as lower limit curve, for the stresses “ »,) was ob- 
tained. 

In connection with the results of the writer’s tests (single circles in Fig. 
8) he assumed the upper limit curve to be a sufficiently accurate bound for 
the limit stress Mm. The limit stress um, of symmetrically situated sec- 
tions is bounded by the pedal curve of this limit curve (Fig. 6). 

The writer repeated the same investigation on the basis of the Mises cri- 
terion of plasticity.4 This leads to limit curves for “4, and Umsg as shown in 
Fig. 8 (reference 4, page 35). The writer’s test results with single and sym- 
metrical oblique sections between holes are given by the single and double 
circles, respectively. Since the rotation of the forces F through an angle y 
(compare Fig. 3) may overstress other sections and since it is too difficult 
to prescribe where the curve for & ms can be used and where not, in refer- 
ence 4 the writer recommended the limit curve for “mp for use for all sec- 
tions. It is an ellipse with the equation 


0243.37? (6) 


It leads to allowable total stresses in the section MM between the holes (Fig. 
5) of (reference 4, Equation (18)) 


Pall (1) 


Um = 
V1 +2.3 sin’ B 


where @,)) is the allowable stress for pure tension, so that the allowable 
force Fis (Fig. 5) (MM)t Um * d) tum. 

The conditions used by the authors for obtaining the upper and lower limit 
loads are essentially identical to those used by the writer in reference 2 and 
lead to the same results. The reasons why the limit loads obtained by the 
authors are not as close together as those in reference 2 is that they did not 
consider the most unfavorable or the most favorable cases. 

The writer does not agree with the authors that the s2/4g specification is 
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rather reasonable, since for a section subjected to pure shear, where g = 0, 
s2/4g would become infinite. Actually their Equation (33), which is based on 
their Equation (31), identical with the writer’s Equation (1), and which 


equates V1 + s2/b2 to 1 + s2/2b2, does not hold if s/b is more than approxi- 


mately unity. If s/b is large, so that V1 + s2/b2 is about equal to s/b, from 
Eq. (33) 


= ktb + 1] = kts 


Actually for large angles § it is more unfavorable to consider the section MM 
between the holes (Fig. 9). If 8 becomes 90°, from Eq. (1) 


(s-d) 


so that the limit stress for 8 = 90° is k, as also follows from Fig. 4. How- 
ever the writer again wants to recommend Equation (7) for the allowable 
stress in the sections MM between the holes (Fig. 5), which, in terms of the 
pitch or stagger s and the gage g, may be written as 


he? all (8) 


s +¢° 


and which he has used in design for about 25 years. 

For more details the writer has to refer to his papers. In his theoretical 
investigations he assumed that the rivets do not keep the holes round, since 
they do not completely fill the holes, as is confirmed by tests done by others 
at that time and later. 


W. G. BRADY® and D. C. DRUCKER, A.M. ASCE.’—The discussion by 
Professor Munse is a very helpful addition to the paper emphasizing as it 
does the limitations of limit analysis. In most, but not all, problems of im- 
portance in civil engineering first yielding or purely elastic action is not real- 
ly of great significance. In most cases the limit load will be the maximum 
load which the engineer would wish to consider and so will base his factor of 
safety upon it. As Professor Munse points out, however, ultimate strength 
may be of more interest in many problems. As to the example of the I-beam 
loaded through the flanges only, it is likely that limit analysis will show an 
efficiency of the order of 75% rather than 100 because work hardening is not 
taken into account. 

Professor Bijlaard’s contributions to the theory of plasticity and its appli- 
cations are well known and it is regretted that no reference was made to his 
papers on the effect of holes and notches. Their presentation in English is 
most welcome. Professor Bijlaard demonstrates the well known fact that an 
experienced engineer can obtain good answers to the most complicated prob- 
lems. It should be noted, however, that the upper or lower bounds he de- 
scribes are not necessarily either. The limit theorems as given in the paper 


6. Cornell Aeronautical Laboratory, Inc., Buffalo, N. Y. 
7. Prof. of Eng., Brown Univ., Providence, R. IL. 
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require that kinematic patterns be employed for upper bounds and that 
equilibrium patterns below yield throughout be employed for lower bounds. 
As reported in the paper the theory and tests both show clearly that the num- 
ber of rows of holes has a significant effect. A simplified discussion in terms 
of what happens in the two row case is not sufficient. Fig. 7 of Professor 
Bijlaard’s discussion would have to be extended to the boundaries of the sheet 
under tension without altering the calculated axial pull. The difficulty in doing 
this, and it is a physical as well as a mathematical one, is precisely what 
causes lower bounds in the paper to be low. The added effort to narrow the 
bounds is very great and not essential in the majority of applications. 

Perhaps the real contribution of limit analysis is not in aiding the engineer 
who knows the practical answers but rather in its value as a tool for the in- 
experienced to produce results which otherwise take a lifetime of learning. 
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Discussion of 
“DESIGNING ALUMINUM ALLOY MEMBERS FOR. 
COMBINED END LOAD AND BENDING” 


By H. N. Hill, E. C. Hartmann, and J. W. Clark 
(Proc. Sep. 300) 


HERBERT A. SAWYER, yr.,} A.M. ASCE.—Certain important refinements 
to present practice in designing beam-columns have been presented in this 
paper, these refinements applying not only to aluminum members, which were 
exclusively considered by the authors, but, as the authors hint and as the 
writer believes, to members of other structural materials as well. 

The authors present evidence and conclude that interaction Eq. 9 gives a 
reasonably accurate, simple means for calculating the strength of aluminum 
alloy beam-columns. The authors’ justification of this equation for failures 
by bending in the plane of the applied bending moment was entirely empirical, 
and therefore they carefully state that additional investigative work is needed 
to determine the limits of this equation’s applicability for conditions differing 
from those of the cited tests. The cited tests were restricted to (1) aluminum 
alloy members with (2) applied bending constant throughout length and with 
(3) applied bending stresses less than half of the ultimate bending modulus 
(fp < fp/2). 

The purpose of this discussion is to consider certain theoretical and ex- 
perimental evidence on the use of this equation outside of these restrictions. 


Theoretical Evidence 


The writer has previously derived? an expression for the strength of beam- 
columns very similar to Eq. 9: 


f 


where {co #2 strength per unit area of a short column. 
2 
= n ae where E; is the tangent modulus of the material at a 
( 
stress of f, (not fc). 

K =a constant which depends on the distribution of the applied 
moments, that is, the distribution of either initial eccentricities 
or initial moments from transverse loads; for triangular distri- 
bution equal to about -0.18, for sinusoidal 0.00, for parabolic 
+0.03, for third-point trapezoidal + 0.05, for quarter-point 
trapezoidal +0.13, and for rectangular about + 0.26. 


1. Assoc. Prof. of Civ. Eng., Univ. of Connecticut, Storrs, Conn. 
2. “Applied Column Theory” by F. R. Shanley, discussion by H. A. Sawyer, 
Transactions ASCE, Vol. 115, 1950, p. 743, Eq. 71. 


606-23 


_| 


This equation was derived from the basic column equation, f = P/A - Mc/I, 
using a generalized secant-formula type of derivation and assigning to “f” a 
sliding, limiting value varying from the short column strength for concentric 
loadings to the modulus of bending (or rupture) for pure bending loadings. Of 
course the main theoretical limitation of Eq. 12 is its inaccuracy above the 
elastic limit; the use of Ey instead of the elastic modulus somewhat reduces 
this inaccuracy. Another limitation of Eq. 12 is that it applies only to failures 
from bending in the plane of the applied moment. However, this is the case 
for which the authors gave only empirical justification for Eq. 9. The folliow- 
ing comparison of Eqs. 9 and 12 is made to reveal certain theoretical justifi- 
cations and limitations for this case of Eq. 9. 

The factor containing K of Eq. 12 is significant in giving what the writer 
believes to be a fairly correct indication of the relative adjustments which 
should be made in the bending term of Eq. 9 to account for distributions of 
applied moment which are not of the rectangular shape considered in this 
paper. For example, to take the most extreme case, if the primary bending 
of a beam-column is caused by a single transverse load at the mid-length (or 
out-of-lineness caused by a single short bend of the axis at the mid-length) 
the primary moments would be triangular in distribution, and according to 
Eq. 12 the bending term of Eq. 9 should be multiplied by the factor: 


| - 0.18 
+026 


cE 


However, the importance of this factor should not be overestimated; it dif- 
fers from unity significantly only when the term {/fck is relatively large, 
but for large values of f, the bending term of the interaction equation is rela- 
tively small and insignificant, and large percentage variations of this small 
term will have little effect on the allowable axial load. Therefore, for most 
design purposes the omission of this variation in Eq. 9 is justified, especially 
since this omission is always on the safe side. 

The remaining differences between Eqs. 9 and 12 are also more apparent 
than real. The term {CE differs from fog only for large values of f,, but for 
such values the bending term of the interaction equation is relatively small. 
This difference tends to make Eq. 9 slightly less conservative than Eq. 12. 
The remaining difference, that between fc and fco, disappears for low values 
of L/r and is unimportant for low values of {-, but for the higher values of 
L/r and f, it becomes more significant, tending to make Eq. 9 more conserva- 
tive than Eq. 12 in this range. 

It may be concluded that with the exception of the range just mentioned Eq. 
12 lends a sort of theoretical justification to Eq. 9 encouraging the expectation 
that Eq. 9 would apply for materials and loadings other than those of the cited 
tests. 

Any designer attempting to use Eq. 9 (or Eq. 12) should note that the Euler- 
Engesser column equation which gives the basic column strength differs from 
the usual practical design formulas in that it neglects the effect of accidental 
eccentricity and heterogeneity. Also, the effect of accidental eccentricity was 
meticulously eliminated in the cited tests. Since such effects cannot be elimi- 
nated in practical construction, they always exist for both principal bending 
axes. The probable maximum accidental eccentricities must be predicted and 
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the resulting stresses included in the term fp. Provision for accidental ec- 
centricities by increasing the factor of safety of the Euler-Engesser formula 
is not satisfactory because eccentricity, and hence the percentage reduction 

in allowable load from eccentricity, tends to increase with increase in column 
length. For example, a standard out-of-lineness tolerance for aluminum alloy 
tubing is one-twelve-hundreth of its length. 


Experimental Evidence 


The surest way to justify extension of the range of Eq. 9 is by comparison 
with test results. Messrs. Tuckerman, Petrenko, and Johnson5 have per- 
formed extensive tests of steel and aluminum alloy tubing under combined 
axial and third-point transverse loads. The writer’s Fig. 1 shows a compari- 
son between these test results, Eq. 9, and Eq. 12 for tests of duralumin tubing. 
Eq. 9 conforms to these test results very well except for the highest values of 
L/r. This exception agrees with the conclusion of the above discussion. The 
writer believes that the disagreement for L/r = 30 may not be significant be- 
cause the modulus of bending tends to increase for low values of L/r and such 
increases were neglected in the report of these tests. Note that Eq. 9 is on 
the side of safety for all values of L/r. 

Chrome-molybdenum steel tubing was used for the remaining Tuckerman 
tests. The writer has shown elsewhere® the remarkable agreement between 
Eq. 12 and these tests. Eq. 9 may also be compared with these test results, 
and such comparison shows a good agreement for values of L/r from 30 to 
70. As would be expected, Eq. 9 becomes conservative for higher values of 
L/r, the curve of Eq. 9 for L/r * 90 practically coinciding with most of the 
experimental curve for L/r * 100, and the curve of Eq. 9 for L/r = 120 prac- 
tically coinciding with most of an extrapolated experimental curve for L/r = 
140. 


Conclusion 


It may be tentatively concluded on the basis of these comparisons with 
theory and the Tuckerman test data that Eq. 9 may be applied with safety to 
certain beam-columns: (1) of some metals other than aluminum; (2) with ap- 
plied bending moments from transverse loads, and (3) with applied bending 
stresses greater than half the bending modulus. Also, Eq. 9 is superior to 
Eq. 12 for general design purposes: it follows somewhat better, according 
to the writer’s check, the aluminum alloy test data of the paper; its variation 
is generally on the safe side; and, most important, it is suitable for failures 
by “lateral” buckling (buckling about a bending axis perpendicular to the bend- 
ing axis for the applied moment). 

Lateral buckling neglected, the chrome-molybdenum-steel example and 
theoretical considerations would seem to indicate that Eq. 12 would prove 


4. “Alcoa Structural Handbook,” Aluminum Co. of America, Pittsburgh, Pa., 
1945, p. 174. 

5. “Strength of Tubing Under Combined Axial and Transverse Loading,” by 
L. B. Tuckerman, S. N. Petrenko, and C. D. Johnson, Technical Note No. 
307, National Advisory Committee for Aeronautics, Washington, D. C., 
June, 1929. 

“Applied Column Theory” by F. R. Shanley, discussion by H. A. Sawyer, 
Transactions ASCE, Vol. 115, 1950, pp. 717, 744, Figs. 15, 23. 
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more precise than Eq. 9 for materials such as mild and certain alloy steels, 
wood, and cast iron (with either sharp knees or somewhat uniform curvature 
in their stress-strain diagrams). However, on the basis of the presently 
available data, Eq. 9 seems to be the more generally applicable and safe, as 
well as being simple and suitable for design. The authors’ Fig. 4 demon- 
strates the inadequacy of the conventional interaction formula for beam- 
columns. Therefore the authors of this paper have made a significant and 
valuable contribution to the art of structural design. 


H. N. HILL,7 M. ASCE, E. C. HARTMANN,® M. ASCE, and J. W. CLARK,? 
J.M. ASCE.—The discussers have made several interesting observations on 
the design of members for combined end load and bending. Mr. McCalley’s 
equation (a) expresses the critical bending moment for the general case in 
which both the cross section of the member and the condition of loading are 
unsymmetrical. Equation (a) represents the same relationship as that de- 
rived by Goodier (reference 9, equation [9]), but Mr. McCalley has succeeded 
in stating the relationship in a manner that is particularly useful for design 
purposes. A special case that might be of interest is the case in which the 
cross section of the member is symmetrical but the loading is not; that is, 
the plane of the applied bending moment does not coincide with one of the 
principal axes. In this case “e” is zero and equation (a) reduces to the follow- 
ing: 


Mer KL Vly GI 


It is worth noting that the value of the end fixity factor “K” in the foregoing 
equations always depends upon the restraint against rotation of the end cross 
sections about the minor principal axis, even though in the case of unsymmet- 
rical bending the minor principal axis is not in the plane of the applied bend- 
ing moment. Thus, the ends of a beam may be “fixed” in the lateral direction 
(normal to the plane of the applied bending moment) yet not fixed in the direc- 
tion that is of interest from the standpoint of determining “K.” This was dem- 
onstrated in a recent paper on the strength of aluminum alloy Z beams sub- 
jected to pure bending moment. 10 

It is gratifying that Professor Sawyer has furnished additional evidence in 
support of the use of equation (9) in the design of aluminum alloy beam- 
columns. His theoretical analysis of equations (12) and (9) supplements the 
discussion of interaction equations and their rational basis in section VII of 
the original paper. 

Prof. Sawyer points out that the Euler-Engesser Column Formula does not 
include the effects of accidental eccentricity or heterogeneity, and he con- 
cludes, “Provision for accidental eccentricities by increasing the factor of 
safety of the Euler-Engesser formula is not satisfactory because eccentricity, 
and hence the percentage reduction in allowable load from eccentricity, tends 
to increase with increase in column length.” It is true that initial crookedness 


“7. Asst. Chief, Eng. Design Div., Aluminum Co. of America, New Kensington, 
Pa. 
8. Chief, Eng. Design Div., Aluminum Co. of America, New Kensington, Pa. 
9. Research Engr., Aluminum Co. of America, New Kensington, Pa. 
10. “Lateral Buckling of Channels and Z-Beams,” by H. N. Hill, Proceedings, 
A.S.C.E., Separate No. 334, November, 1953. 
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tends to increase with increasing length of column; however, it does not follow 
that the factor of safety resulting from the use of a formula for ideal columns 
decreases for the longer columns. As a matter of fact, if eccentricity is di- 
rectly proportional to length, the minimum factor of safety given by a formula 
for ideal columns occurs in the general vicinity of the juncture between the 
elastic and inelastic portions of the column curve. 

The above statement is illustrated in Figs. 9 and 10. The solid curve in 
Fig. 9 is the Euler-tangent modulus column curve for aluminum alloy 61S-T6 
rectangular tubes tested in the investigation described in reference 16. The 
lower, dashed curve in Fig. 9 represents the column strength calculated from 
equation (9) for the same specimens loaded as pin-ended columns with a ratio 
of eccentricity, “A”, to length, “L,” of 1/1200. None of the column specimens 
described in reference 16 had a *4/L” ratio of exactly 1/1200. However, test 
values of column strength corresponding to this ratio of eccentricity to length 
were found by interpolation from the test data. These interpolated test values 
are also plotted on Fig. 9. 

Fig. 10 shows the ratio of column strength for a “ 4/L” of 1/1200 to the 
strength of an ideal column determined from the curves and test points in 
Fig. 9. Both the calculated curve and the test points in Fig. 10 illustrate the 
fact that the reduction in column strength caused by eccentricity does not 
necessarily increase with increasing length, even though the eccentricity it- 
self does. Therefore, it does not appear unreasonable to design columns on 
the basis of formulas for ideal members, taking account of the effects of ac- 
cidental eccentricity in the factor of safety. 


April 16, 1954 


t 
H 
606-28 


EXPERIMENTAL VALUES 
OF fe FOR RECTANGULAR 
TUBES WITH = 


(INTERPOLATED FROM 
DATA IN REFERENCE 16) 


WALL 
RECTANGULAR 
TUBE, 61S-T6 
ALLOY 


FOR 
TUBES WITH t * i200 


WwW 
a 
4 
x 
x 
z 
= 
” 
” 
W 
a 
= 
oO 


SLENDERNESS RATIO, 


FIG. 9 - STRENGTH OF COLUMNS WITH —& 


60 
40 
30 L 
fe FROM EQUATION (9) \ | 
20} 
| re) 
@) 20 40 60 80 100 120 
1200 
606-29 


f¢/te FROM EQUATION (9) FOR 61S-T6 
RECTANGULAR TUBES WITH 4 s 7355 


O EXPERIMENTAL VALUES OF 
FOR 61S-T6 RECTANGULAR TUBES 


SLENDERNESS RATIO + 


FIG. 1O- REDUCTION IN COLUMN STRENGTH WHEN £ 


* 7200 


0 20 40 60 80 100 120 
606-30 


Discussion of 
“PLASTIC DEFORMATION OF WIDE-FLANGE BEAM-COLUMNS’” 


By Robert L. Ketter, Edmund L. Kaminsky, and Lynn S. Beedle 
(Proc. Sep. 330) 


ROBERT L. KETTER! J.M., and EDMUND L. KAMINSKY,” J.M., and LYNNS. 


BEEDLE,” A.M. ASCE.—The discussion by Professor Winter focuses attention 
on several points of importance. First of all, each of the members tested as 
part of the experimental study was loaded with a combination of axial thrust 
and end moments applied about the stronger principal axis of the section. 
Moreover knife edges were located at each end of the specimens in a manner 
which allowed rotation about an axis perpendicular to the web but restrained 

it in the direction of the web. This was done since it was considered that 
strong-axis bending is that condition most frequently encountered in practice. 

From the previous Lehigh study cited by Professor Winter? it is possible 
to compare the predicted strength of columns with the results of pin-end tests 
in which the possibility of twist was eliminated. In Figure 42 is shown the 
correlation between the predicted strength curve (Figure 38, ec/r2 * 1.0) and 
a set of weak axis tests (315.7). The collapse solution neglecting residual 
stress has been used even though the members were tested in the as-delivered 
condition. This was done since the pure axial test of the same section indi- 
cated that a relatively small amount of cooling residual stresses was present 
in the specimen. Agreement between the predicted and observed values is 
reasonable. 

A few tests have now been made on columns of 8WF31 shape bent about the 
weak axis. The results will be included in a future report; they are in excel- 
lent agreement with the theory. 

The second point raised was that of the generality of the assumed cooling 
residual stress pattern (Figure 7). Although the actual magnitude and pattern 
of residual stress present in any given cross-section is influenced by its geo- 
metric proportions, stresses observed in the flanges of column-type sections 
actually do not vary as much as might be expected. 

In Figure 43 is plotted the distribution of residual stresses in four different 
WF columns. They indicate that the assumption made in the paper was not un- 
reasonable. All of these shapes have compressive residual stresses at the 
flance edges. More recent measurements made on still other WF column 
shapes (12WF50, 12WF65, 6WF15.5, 14WF43) show a similar result. 

Since the investigation is concerned with columns and not with beams, then 
shapes with flanges that are narrow in comparison with beam depth are not of 
prime interest. Further the 8WF31 shape is rolled in more tonnage than most 


1. Research Asst., Lehigh Univ., Fritz Eng. Lab., Bethlehem, Pa. 

2. Dept. of Civ. Eng., Syracuse Univ., Syracuse, N. Y. 

3. Asst. to the Director, Lehigh Univ., Fritz Eng. Lab., Bethlehem, Pa. 

4. Johnston, Bruce and Cheney, Lloyd, “STEEL COLUMNS OF ROLLED WIDE 
FLANGE SECTION”, American Institute of Steel Construction, Progress 
Report #2, November, 1942. 
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other columns. Finally, since the flanges provide the majority of the bending 
stiffness, the magnitude and pattern of residual stress present in the web of 
the specimen is of relatively minor importance. 

Professor Winter states, “The quoted test results seem to indicate that 
this feature, influence of twist, may prove tobe of far greater consequence in 
reducing the strength of such members than the influence of residual stresses 

...”, such a separation of the influence of twist from that of residual stress 
may be misleading. Torsional stiffness is made up of essentially two parts: 
a St. Venant type of torsional stiffness and a cross-bending type. Discounting 
the first factor for the purposes of this discussion, (its influence will decrease 
as length decreases) the second factor, cross-bending stiffness, is markedly 
affected by residual stress. If, for example, a flange had yeilded due to re- 
siduals and applied stress to a distance equal to 1/4 of the width of the 
flange from each side, the cross-bending stiffness of that flange containing 
residuals would be 1/8 as large as for the same one without such initial stress 
imperfections. True, the inelastic lateral-torsional problem has not been 
solved and a solution is urgently needed, but this solution must include the in- 
fluence of residuals if it is to solve the problem that is present in practice. 

The final point concerns the correlation with the restrained column tests 
carried out at Cornell (see Figure 41). As was discovered independently both 
by the authors and the discusser, the implied application of the theory to this 
particular problem does not lead to the correct solution. Instead, the proper 
procedure is as follows: 


1. Given the spring constant, 8, assume effective length for various P/Py 
values. 
2. Determine 7 from Figure 44 knowing py and cL/r. 


(Note: 1 is the secant to the M/My - § / $y curve at the point corre- 
sponding to collapse of the member in question. This information was 
available from the calculations leading up to Figures 37 and 38.) 
Check “c” by using the Cornell 8o/n EK versus “c” curves.® 

. Repeat steps 1., 2., and 3. until the assumed “c” values agree with those 
calculated. 
Using the correct “7” values corresponding to each assumed P/Py, 
determine the relation between € and e from Eq. 8 of the Cornell Report 
(ASCE Sep. 292). 

. Determine (P/ Py) and (€/r) for the corresponding (cL/r) values from 
Figure 38 of this report. 

7. Obtain e/r versus P/Py by multiplying the results of steps 5. and 6. 


The corrected curve is shown in Figure 45. It is in reasonable agreement 
both with the test results and with the Cornell analytical solution. 


5. This procedure is a modification of one developed by Messrs.: Bijlaard, 
Fisher, and Winter and has been presented in more detail in the authors’ 
discussion of ASCE Separate 292 (Vol. 79). 

6. Winter, G., Hsu, P. T., Koo, B., and Loh, M. H., BUCKLING OF TRUSSES 
AND RIGID FRAMES, Cornell University, Engineering Experiment Station, 
Bulletin 36, April 1948. 
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CORRECTION TO DISCUSSION OF THE APPLICATION OF 
HEAVISIDE’S STEP-FUNCTION TO BEAM PROBLEMS 


(Proc. Sep. 202) 


On page 333-17, in the fifth line from the bottom, the coefficient 0; should 
be changed to Cj. On page 333-18, in Eq VI, the factor (x - 2) should be 
changed to (x - a). To the line following Eq VI, should be added the term: 
+bg k sen k (x - a). Also on page 333-18, in line 25, the last term within 
brackets should be ag instead of a. 
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Discussion of 
“DEFLECTIONS OF A CIRCULAR BEAM OUT OF ITS INITIAL PLANE” 


by Enrico Volterra 
(Proc. Sep. 375) 


I. OrsterBLom,? M. ASCE.—In the English language, the subject of the 
paper was ably treated as early as 1914 by A. H. Gibson and E. G. Ritchie (22). 
Two distinguished authorities, Messrs. Pippard and Barrow, made extensive 
studies of and published papers and books on the arc bow girder (10) (11). 
One could ask for no more complete treatment of the theory of bow girder 
deflections and of the application of this theory. A book (23) which is a 
résumé of many years work has been published by Messrs. Pippard and 
Barrow. The German references are also extensive, but none that specifically 
refer to the subject or have usefulness in design was published earlier than the 
English works. Thus, one questions whether or not the paper advances new 
theories in the design or construction of civil engineering work, especially as 
the selected solutions do not include a graph for the almost universal case, 
2 y = 180°. Arc bow girders are usually supported at right angles to the 
structural frame and, when they are close to the frame (these are usually only 
ornamental in purpose), the designing engineer need not—and cannot afford 
to—spend time on extensive calculations. The author suggests (Fig. 1) five 
equidistant cuts of the circle, without explanation in the text, and develops 
three of them into helpful graphs. However, if wall wedges of 120° and 45° 
were common, the arc bow girders would be placed not at the point of the wedge 
but somewhere on the face of the wall so that 2 y would equal 180°. 

Solutions to more comprehensive, varied, and real problems than those 
of this paper may be found in papers and books previously published. The 
author has referred to some that seem superior to his own presentation from 
the standpoint of the progress indicated and the utility shown. A bulletin 
(24) by Mr. Moorman provides an objective understanding of past work 
in this field, an experimental verification of the theory, and (best of all) a 
clarification of what the engineer may expect and should demand from the 
mathematician. 

The paper fails to elaborate on the relation of stress to deflection ; yet, because 
of the combination of longitudinal and rotational deformations, this subject is 
important, particularly in reinforced concrete design. The practicing engineer 
would stand defeated at this point even if the author’s graphs would otherwise 
suit his needs. The paper would have been much more important if the author 
had attempted such an interesting problem. 

In Appendix II, one of the books (19) seems to have no reference to the 
subject; two of them (3) (20) refer to mechanical engineering problems; and a 
fourth (1) is so broad and fundamental that it is the basis for all elastic problems. 
The writer does not consider these works proper for inclusion in the Biblio- 
graphy on Curved Beams. 


Enrico Votrerra,* M. ASCE.—The writer thanks Mr. Oesterblom for 
contributing to the discussion of the paper. 


3 Rutherfordton, N. C. 
* Prof. of Mechanics, Rensselaer Polytechnic Inst., Troy, N.Y. 
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The purpose of the paper was the study of deflections of circular beams in 
three particular and distinct cases: (1) Circular rings rigidly supported at 
symmetrical points and subjected to uniformly distributed forces; (2) circular- 
arc bow girders built in at both ends and subjected to uniformly distributed 


forces; and (3) circular-arc bow girders built in at both ends and subjected to 
concentrated forces. 


Fig. 1 refers to case 1 of the circular ring and not to case 2 of the circular- 
arc bow girder (Fig. 2(a)) as erroneously interpreted by the discusser. 
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Originally, the writer prepared many more graphs than those that were 
finally published; among the unpublished graphs were those for the case of 
circular-are bow girders with a center angle 2-y = 180° to which Mr. Oester- 
blom refers. Graphs for this case are presented in Figs. 12 and 13. 

The writer is aware of the works previously published on the subject, 
particularly those of English authors. He presented a bibliography which he 
thought was almost complete, containing all the works published on curved 
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beams from the 1843 papers of Saint-Venant to the papers published on the 
subject in the United States in 1950. However, since the presentation of the 
paper in the spring of 1952, other works have appeared on the subject, among 
them the excellent textbook of Mr. Pippard (25). In chapter nine of this 
book the problems of bow girders and transversely loaded frames are dis- 
cussed (25a). Mr. Pippard’s approach to the problems is as follows. He 
expresses the total strain energy due to bending and twisting, neglecting 
the small contribution to the total energy due to shear. By applying the 
principle of minimum strain energy, the three redundant unknowns—the 
bending and twisting moments and the shearing force in the midsection of the 
bow girder—can be computed. In order to simplify the numerical computa- 
tions involved in solving the three simultaneous equations, the method of 
superposition is applied by dividing the loading into a symmetrical and skew- 
symmetrical system. These two cases are dealt with separately, and the 
results are summed algebraically to obtain the final required solutions. 

By this method, instead of solving three simultaneous equations, the 
solution of a first equation gives the unknown bending moment whereas the 
solution of a pair of simultaneous equations yields the unknown twisting 
moment and shearing force, which can therefore be obtained explicitly. 

Mr. Pippard discusses many cases of arc girders of engineering interest, 
among which are the cases of a circular-arc bow girder with a concentrated load 
and a circular-arc bow girder with a distributed load. 

However, in the more general cases of noncircular-are bow girders, Mr. 
Pippard finds that the resulting expressions 

“*** cannot usually be integrated by exact mathematical process but 

must be calculated approximately by graphical or arithmetical methods.” 

In the paper, the writer used a different and, he believes, a more general 
method of approach to determine bending and twisting moments, M, and M,, 
deflection, v, and twisting angle, 8. Four differential equations, connecting 
the four quantities, M,, M,, v, 8, the initial curvature of the central line of 
the beam, R, and its inertial and torsional characteristics, E J and K, respect- 
ively, are used. Two of these equations are the well-known Saint-Venant 
equations (Eqs. 1) whereas the two subsidiary equations of equilibrium (Eqs. 7) 
are obtained by minimizing the potential energy of the system. In the 
equation for the potential energy only the terms expressing bending and 
twisting moments appear because, as Mr. Pippard indicated, the energy caused 
by shear can be neglected. By solving these four differential equations under 
different load and boundary conditions, the results for the three special cases 
mentioned were obtained in closed form. 

The same differential equations under appropriate load conditions and 
boundary conditions were used in earlier papers by the writer for solving 
problems of circular rings resting on elastic foundations and submitted to 
symmetric or antisymmetric concentrated loads (2) (21) (26). 

Moreover, by using the same approach it is possible to solve in closed form 
particular cases of the more general problem of the deflection of a noncircular 
bow girder. In fact, it can also be shown that, in the case of a noncircular 
bow girder built-in at both ends and subjected to a uniformly distributed load, 


606-42 


the results can be expressed in closed form if the central line of the curved 
beam is a cycloid, a catenary, or a parabola. The writer is preparing a paper 
(1954) on the subject of deflections of noncircular bow girders in which the 
results for these special cases will be given in closed form. 

Mr. Oesterblom’s criticism that the paper fails to elaborate on the relation 
of stress to deflection is certainly not correct. In fact—as the formulas pre- 
sented in the paper give in closed form the relationships existing among bending 
and twisting moments, deflections and angles of twist, initial curvatures, 
inertial and torsional characteristics of the beam, and external loading—it is 
possible in each specific case to derive from these formulas the relationships 
between stresses and deflections. 

The textbook (19) to which the discusser refers as having no reference to 
the subject contains an entire chapter (Chapter XXI) dealing with small 
deformation of naturally curved rods (19a). 

The two references (3) and (20), which according to the discusser apply to 
mechanical engineering problems and should not have been quoted, are Mr. 
Resal’s fundamental work on bending and twisting of a curved beam and the 
textbook of 8. Timoshenko, many pages (of the second volume) of which 
deal with the problem of bending of a curved beam out of its plane of initial 
curvature and give useful references on the entire subject. 

Finally, the fourth work, which Mr. Oesterblom does not consider proper 
for inclusion in the bibliography on curved beams, is Saint-Venant’s origina] 
papers (1). In these papers the theory of small deformations of initially 
curved beams was first developed. 
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Discussion of 
“AN EXPERIMENTAL STUDY OF BUBBLES MOVING IN LIQUIDS” 


By W. L. Haberman and R. kK. Morton 
(Proc. Sep. 387) 


J. A. ROBERSON, ! C. L. BARKER,” A.M. ASCE, and P. F. RUFF,? J.M. 
ASCE.—Messers. Haberman and Morton are to be complimented for their 
work concerning the movement of bubbles through liquids. Certainly their 
paper will be a definite contribution to the literature on this subject. 

At the State College of Washington a study has been made by the writers 
of the movement of bubbles in a turbulent stream. It is interesting to note 
that close agreement exists between the results reported by the authors and 
those observed by the writers. 

The problem of determining the rise velocity of bubbles, although some- 
what academic to many becomes of great importance, e.g. in pipe lines where 
air is mechanically entrained. In lines of this type small bubbles tend to coa- 
lesce; and eventually form large air masses or “pockets.” The removal of 
these air masses by pipe vents produces pressures that can be destructive. 
The problem of locating vents in a line for the removal of bubbles before they 
can form air masses, involved a study of bubble rise in a turbulent stream. 

If U is the rise velocity, V the velocity of the water, and D the pipe diame- 
ter, the travel, L, of a bubble in rising a distance D is: 

._D 
L* v U 
In the study of vent location in pipe lines, L, in the above formula, is the dis- 
tance the vent should be placed downstream from the point of entrance of the 
air bubbles into the pipe line, and U is the minimum rise velocity of the bub- 
bles. 

In the studies conducted by the writers, bubbles 1/4 inch in diameter or 
greater were studied. Davies’ found the minimum rise velocity of bubbles 
of this size range, in a stationary liquid to be 0.6 feet per second. However, 
the problem with which the writers were concerned involved a liquid in a state 
of turbulence. An experimental apparatus was constructed of 5 1/2 inch lucite 
tubing to study the rise velocity of bubbles in a turbulent stream. The rise 
velocity, U, in a turbulent stream (Reynolds Number * 50,000 to 300,000) as 
determined by the writers agrees closely with the results obtained by Davies 


1. Asst. Hydr. Engr., Div. of Industrial Research, and Instructor, Dept. of 
Civ. Eng., The State College of Washington, Pullman, Wash. 

2. Prof. of Hydr. Eng., The State College of Washington, Pullman, Wash. 

3. Junior Hydr. Engr., Div. of Industrial Research, and Instructor, Dept. of 
Civ. Eng., The State College of Washington, Pullman, Wash. 

4. “The Mechanics of Large Bubbles Rising Through Extended Liquids and 
Through Liquids in Tubes”, by R. M. Davies and G. I. Taylor, Proceedings, 
Royal Society of London, Vol. 200, Series A, 1950. 
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and the authors. Precise comparison was not practical because of the varia- 
tions above and below the mean rise velocity as a result of turbulence. 

It was of interest to the writers to note that the data of Figure 2, of the 
authors’ paper, plotted a smooth curve. The results presented by Davies, as 
did the writers, gave a spread of points, especially for the bubbles of relative- 
ly large diameter. 
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CORRECTION TO DISCUSSION OF APPLICATIONS OF THE 
RELAXATION TECHNIQUE IN FLUID MECHANICS 


(Proc. Sep. 223) 


The name of Mladen Boreli was misspelled on the inside front cover and 
on page 396-3. 


CORRECTION TO DISCUSSION OF MECHANICS OF MANIFOLD FLOW 


(Proc. Sep. 258) 


On page 396-22, in the first term of the right-hand member of Eq. 4, the 
Squared quantity should be summed to i for each value of i before summing 
the entire system to n-1. 
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